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PROGRAM SUMMARY 



Title of program: SU3 
Catalogue number: 



Program obtainable from: CPC Program Library, Queen's University of Belfast, 
N. Ireland (see application form in this issue) 

Licensing provisions: Persons requesting the program must sign the standard 
CPC non-profit use license (see license agreement in every issue). 

Computers for which the program is designed and tested: DEC VAX 4000/90, 
DEC Alpha 3000/700, Dell 486, Dell Pentium, HP715/33 

Operating systems under which the program has been tested: VMS v6.2, MS- 
D0S6.X, HP UNIX 9.0 

Programming language used: C Language 

Memory required to execute: Variable; on the VAX it requires 934 k words for 
3 (g) 3, 998 k words for 27 27 

Disk space used for output: Variable; 2 k words for 3 (g) 3, 182 k words for 27 
27 



No. of bits in a word: 32 
No. of processors used: 1 



Has the code been vectorized? No 



No. of lines in distributed program: 2601 



Keywords: SU{3), Clebsch-Gordan coefficients, Wigner coefficients, vector cou- 
phng coefficients, isoscalar factors 

Classification: 4.2 Computational methods of algebras and groups 
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Nature of physical problem: Calculations using models based on SU{2t) symme- 
try often require the Clebsch-Gordan coefficients which give the amplitudes for 
the expansion of direct products of irreducible representations. 

Method of solution: Recursion relations [1] are used to construct tables of the 
isoscalar factors for the product of two irreducible representations. The SU (2) 
Clebsch-Gordan coefficients [2] [3] are then used with the isoscalar factors to 
produce the SU (3) Clebsch-Gordan coefficients. 

Restrictions on the complexity of the problem: The program is limited by com- 
puter memory and maximum size of the integer variables. On the DEC VAX 
with 32-bit integers, the largest product that has been successfully computed is 
27 ® 27. 

Typical running time: Running time varies according to the sizes of pi, gi, p2, 
and q2 in the product (pi, qi) ® (pa, 92)- On the VAX for 3 3 the CPU time 
is 1.2 sec; for 27 ® 27 the CPU time is 13 min 11 sec. On the Alpha, the latter 
time is only 1 min 13 sec. 

Unusual features of the program: The program has these unique features: it 
determines the representations in the Clebsch-Gordan series; it uses recursion 
relations to efficiently compute the isoscalar factors; it presents exact results; it 
provides results as I^TgX [4] formatted tables. 

References 

[1] H. T. WiUiams, Symmetry Properties of SU3 Vector Coupling Coefficients, 
WLUPY-9-93. 

[2] E. U. Condon and G. H. Shortley, The Theory of Atomic Spectra (Cam- 
bridge University Press, London, 1957). 

[3] E. P. Wigner, Group Theory and its Application to the Quantum Mechan- 
ics of Atomic Spectra, trans, by J. J. Griffin (Academic Press, New York, 
1959). 

[4] L. Lamport, I^TgX: A Document Preparation System (Addison- Wesley, 
Menlo Park CA, 1986). 



2 



LONG WRITE-UP 



1 Introduction 

In calculations in nuclear and particle physics which assume SU{3) symmetry, 
it is often useful to have values for relevant Clebsch-Gordan coefficients (CGC) 
0. Some tables of limited extent have been published in the past f!]] |H [0 Pf 



]. Also, programs have been written which generate these coefficients [T0| ||TT 



1^ [|T^. This new program is meant to replace them. It has some advantages 
that make it faster and easier to use. Comparative speed is achieved by use of 
recursion relations to generate the isoscalar factors (ISF), from which the CGC 
are constructed. This method requires less memory and allows larger tables to 
be generated. All calculations are done exactly, using integer variables. Output 
is in the form of li^TgX [jl4| formatted tables. The code is written in ANSI 
standard C, so it should run on most conventional platforms. 



2 Preliminaries 

The Clebsch-Gordan coefficients are the amplitudes for the projection of the 
product of two irreducible representations (ri, of SU{3) onto the separate 
irreducible representations (R„) found in the Clebsch-Gordan series: 

ri®r2 = 5]R„. (1) 

They can be defined by 

|riai) |r2a2) = (R-i A|riair2a2) \^iA) , (2) 

where Ai, ai, and a2 denote the quantum numbers specifying the particular 
states within the representations. Particle physicists are accustomed to the 
quantum numbers hypercharge y, isospin i, and third component of isospin i^] 
herein are also used the "projection" quantum numbers k, I, m, which are related 

*The Clebsch-Gordan coefficients are also called vector coupling coefficients or Wigner 
coefficients. Readers interested in a good theoretical background are directed to 0] 1^ |^ j|] . 
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to y, i, and iz in a representation (denoted by the representation variables (p, q)) 
by 

k = |(p + 2g) + iy + z, 

I = |(p + 2g) + iy-^, (3) 
m = lip + 2q) + ly + iz, 

The highest-weight state (shw) in a given representation is defined to be that 
which has the largest iz] thus the state with 

1—0 

"■shw 

The (outer) degeneracy of an irreducible representation in the Clebsch- 
Gordan series (Equation |1|) can be determined algebraically from the represen- 
tation variables [|15|. A representation appears in the series if the degeneracy 
is greater or equal to one, and only in such cases does the program attempt to 
calculate its coupling coefficients. If we represent the irreducible representations 
as 

ri = (PuQi), 

r2 = (P2,g2), (5) 
R = {p,q); 



and define 

then if both 
1. a G 



a = I [{Pi - Qi) + {P2 - q2) -{p- q)] , 
b = |[(pi + 2gi) + (p2 + 2g2)-(p + 2g)]; 



(6) 



Z (integers) (thus 6 G Z), and 

2. < 6 < min(gi + q2,pi + qi,P2 + ^2), -min(gi,g2) < a < min(pi,p2), 
< a + 6 < min(pi + qi,p2 + q2,Pi +P2), 

the representation appears in the series and has degeneracy 

d{p,q:pi,qi;p2,q2) = 1 

+ mm{q2,pi + qi,b,pi - a) 

— max(0, b — qi,b — P2, —a, b — a — qi,a + b — ^2)- 

(7) 
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The isoscalar factors F are defined by 

(R F / J^lri t/i ii r2 ^2 «2 «2^) = -F(R, F, / : ri, yi, ii; r2, ?/2, ^2) 

/T T I- ■ ■ ■ \ 

X (/ Iz\^l ilz i2 l2z) , 

wliere the last bracket is an SU{2) CGC, which in the Condon-Shortley phase 
convention is [|T^ [|1^ 

(/ Iz\ii iiz ^2 kz) = 
S{Iz - Hz - ■i2z) 



X 



(J + 21 - t2)\il -ti + ^2y■i^l + ^2 - /)!(/ + Iz)K2I + 1) 
\il + ii+i2 + iy.{ii - iizV-in + Hz)Ki2 - ^2.)!(^2 + i2zV. 

^ ^ {I - ii + i2 - n)\{I + Iz - n)\n\{ii -i2-Iz + n)V 
where n runs from 

nio„ = max{0,iiz-ii,Iz-ii-i2) (10) 

to 

rihigh = min(/ + 12 + hz, I - h + i2,I + h)- (H) 

Note that 

(/ IzVi iiz i2 i2z) = i-iy-'^-'^ (J 1,1^2 i2z ii Hz) . (12) 

2.1 Conventions 

The internal phase convention, fixing totally the relative phases between states 
within a representation, is chosen to be that of de Swart The irreducible 
representations of SU{?>) can be thought of as consisting of SU{2) multiplets 
[isomultiplets), each at a specific hypercharge. This internal phase convention 
corresponds to adopting the Condon-Shortley [0] phase convention for the T- 
and IZ-spin operators. (In the fiavor-S'f/(3) notation, interchanges u and d 
quarks, while interchanges u and s.) 

The overall phase of the representations in the product is also chosen fol- 
lowing de Swart [|I[| . Consider the state of highest weight in the product repre- 
sentation, |R shw); the state of highest weight in the first factor, |ri shwi); and 
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the state in the second factor r2 = {p2, (I2) with the highest isospin that couples 
to shw and shwi, |r2 /c2max(shw, shwi) Z2min(shw, shwi)). The outer phase con- 
ventions requires that the Clebsch-Gordan coefficient for couphng these states 
be real and positive: 

(R shw|ri shwi r2 /i;2max(shw, shwi) Z2min(shw, shwi)) > 0. (13) 

The Condon-Shortley convention for the SU2 CGC assures that the correspond- 
ing isoscalar factor is also real and positive. The internal and external phase 
conventions guarantee that all isoscalar factors and Clebsch-Gordan coefficients 
are real. 

In couplings where the degeneracy d > 1 there are d distinct sets of CGC, 
and one must chose by convention a technique for resolving the outer degeneracy. 
Several distinct techniques for this resolution are described in the literature 
18| ||19|| . The technique chosen here considerably simplifies the determination 



of isoscalar factors through use of recursion relations, and produces ISF and 
CGC which share the symmetries under interchange of representations (Racah 
symmetries), which are familiar from the SU2 case. The technique relies on the 
fact that there are nonvanishing ISFs of the form 

F(R, shw : ri, shwi; r2, ^2, ^2) (14) 

for at least as many distinct values of k2 as the degeneracy d. If one chooses ISFs 
for the d—1 highest values of ^2 to be zero, the recursion relations, normalization, 
and the above-mentioned sign conventions uniquely determine a full set of ISFs. 
A second set can be determined by choosing the ISFs for the d — 2 highest values 
of /c2 to be zero, and insisting that this second set be orthogonal to the ffist. 
Subsequent sets of ISFs are formed likewise by successively forcing fewer of these 
ISFs to zero, and enforcing the condition of orthogonality to those previously 
determined. 

This technique of resolution of the outer degeneracies insures that the cou- 
pling coefficients in all cases change by at most a sign under the interchange of 



ri and r2, in contrast to the "canonical" resolution scheme 19 . 
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2.2 Symmetry Phases 



With the conventions of the previous section, all the ISFs and CGCs for SU (3) 
are fully determined. The resulting symmetries are useful in reducing the num- 
ber of values which must necessarily be tabulated. Three phases are involved in 
the symmetries and are defined in 

In the case of interchange of the two factor representations, 

F(R,y,J:r2,i/2,^2;ri,i/i,2i) = {-1)'-'^-'^^^ ■ r^r^) ^^^^ 

F{R,Y,I : ri,yi,ii;r2,?/2,«2), 

where the factor (— l)'^^*i~*2 comes from Equation |12[ The phase ^i(R : ri; r2) 
depends only on the identity of the representations and on our phase conven- 
tions. Consider the reversed product r2 (S> ri. Suppose that the highest-isospin 
multiplet in ri that couples to the state of highest weight in r2 and the shw of 
R with nonzero ISF has quantum numbers y[ and i[. Then 

ei(R : ri; r2) = (-l)^h-i'i-'2h x sign [F(R, shw : r,,y[, z[- r^, shw2)] , (16) 

where sign(a;) = a;/|a;|. 

Should the second factor and the product representations be interchanged, 
a multiplicative phase and a term related to the ratio of the representations' 
dimensions is generated: 

F{r2,-y2,i2 ■■ ri,yi,ii;R,-Y,I) = 



X, 



(P2 + i)(g2 + i)(P2 + g2 + 2)(2/ + i) 



-\ (p+l)(g+l)(p + g + 2)(2^2 + l) 
xF(R, y, / : ri, yi, ii, r2, ?/2, ^2)- 

If all the representations are conjugated, then a phase ^3 enters: 

F(R,r,/:ri,yi,^i;f2,|/2,22) = (-l)^-*^-'^e3(R : ri; r2) 

xF(R,-F,/ : ri,-yi,ii; r2,-?/2,«2) 



(17) 



Consider the isomultiplets described in Section 2.1. Then from Equation |T2 
find that 



we 



6(R:ri;r2) = (-l)^h-hh-^.. (19) 
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3 Calculational Method 



3.1 Recursion Relations for Isoscalar Factors 

Recursion relations among the CGC for SU (3) can be generated by applying one 
of the SU (3) ladder operators to a pair of SU (3) single particle states coupled 
(via CGCs) to a state of good SU (3) quantum numbers. Using the linearity of 
the ladder operators (e.g. = Ki^ + K2+, where is the ladder operator 
for the coupled state, and and K2+ are the corresponding operators for 
the 1 and 2 states), and orthogonality of states with different SU{3) quantum 
numbers, one can find four-term recursion relations for the CGCs. Making use 
of explicit analytic forms for the SU2 Clebsch-Gordan coefficients these can be 
turned into four-term recursion relations for the SU{3) isoscalar factors. Two 
of the ladder operators (/+ and J_) generate recursions only within the SU2 
variables and produce only identities for the ISF. From the remaining four ladder 
operators expressions fully sufficient to generate all the SU{3) ISFs have been 
derived. 

The simplest language for presenting the recursion relations is to use {p, q) 
to indicate an irreducible representation, and {k,l,m) to specify a state within 
the representation. The scheme involves, as a first step, using one of the following 
relations to establish values for all the isoscalar factors representing coupling to 
the state of highest weight (i.e. k = m = p + q,l = 0). A simplified notation 
which omits the p's and g's, and which uses the symbol shw to represent k = p+q, 
I = will be adopted. Thus, for example, 

F{p,q,k = p + qj = : pi, gi, fci, /i; P2, ^2, ^2, ^2) = -F(shw : kiji, k2,l2)- (20) 

The two recursion relations involving the states of highest weight are 

= aiF(shw : /ci - 1, /i; A;2, /2) + a2-^(shw : fci, /i; A;2 - 1, ^2) ^^^^ 
+a3F(shw : ki, k - 1; k2, h) + a4F(shw : fci, ^2, h - !)• 

where 

ai = {ki + l)(A;i - qi){pi + qi - ki + l){p + q + ki - h + k2 - k + 3) 

2{p + q + ki-h-k2 + l2 + 1) 
{k,-h){k^-h + l) 



\ 
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02 = ai(l ^ 2), 



03 = — \/^i(gi-/i + l)(pi + gi-/i + 2) 


(22) 


X \l {-P - g + fci - /i + A;2 - ^2 + 1) 






2{p + q- ki + /i + k2 -12 + I) 




{ki-h + l){ki-h + 2) ' 





04 
and 



-03(1 ^ 2); 

-- 6iF(shw : /ci + 1, /i; k2, h) + b2F{shw : ki, k; /cg + 1, k) 
+63F(shw : fci, /i + 1; /c2, /2) + 64F(shw : /ci, /i; A;2, ^2 + 1)- 



(23) 



where 



h = ^(/ci + 2)(A;i - gi + l)(pi + - ki){-p - q + ki - h + k2 - k + 1) 



2{p + q-ki + li + k2-l2 + l) 
'\ {k,-h + l){k,-h + 2) ' 

&2 = -61(1^2), 



(/i + l)(gi - /i)(pi + qi-h + l){p + q + ki-h + k2-l2 + ^) (24) 



2(p + g + /c2-/2-fc2 + ^2 + l) 
^ {ki-h){ki-h + l) 



h = 63(1^2). 

Once the states of highest weight have been determined, two other relations 
are sufficient to step from these to any non-shw state. They are 



F{k, I : ki, h; k2, h) = ciF{k + 1,1 — I: ki, h; k2, h) 
+C2F{k,l- 1 : ki,li - l\k2M) 
+C3F(k,l- 1 : ki,li;k2 - 1,^2) 
+C4F{k, I — 1 : ki, h; k2, ^2 - 1), 



(25) 



where 



Ci — a 



{k + 2){k - q + l){p + q - k){ki - h + k2 - h - k + I) 
{k-l + 2y{k -l + k^-h + k2-l2 + ^) ' 
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2a, 



/i(gi-/i + l)(pi + gi-/i + 2) 
^ (A;i - Zi + 2)(A; - Z + All - Zi + A;2 - ^2 + 4) ' 



{kr -h + l) 



\\{k-l-ki + h + k2-l2 + 2) 



-a 



(fca + l)(fc2 - g2)(fci -h + k2-l2- k + l) 

{k2-l2){k2-l2 + l) 



iP2 + q2- k2 + l) 



\\{k-l-ki + h + k2-l2 + 2y 



a 



\ 



l2{q2 -I2 + l){k -l + ki-h- k2 + l2 + 2) 

{k2-l2 + l)ik2-l2 + '2) 



(P2 + 92-^2 + 2) 



i - 

\ {k - I + ki - h + k2 - I2 + 4.)' 



{k-l + 2f{k -l-ki + h + k2-l2 + 2) 
\l{q-l + l){p + q-l + 2){k-l + ki-h-k2 + l2 + 2)' 

F{k, : fci, /i; k2, h) = diF{k + 1, : fci + 1, h; k2, h) 
+rf2F(/c+ 1,0 : /ci,/i;/c2 + l,/2) 
+d3i^(fc + l,0: fci,/i;fc2,^2 + 1), 



di 



2/3. 



(A;i + 2)(A;i-?i + l) 



^ (A;i-Zi + 2)(A;-Z + A;i-Zi + A;2-Z2 + 4)' 



^ (A; - / + /ci - /i - /C2 + /2 + 2) 



do 



/3. 



{k2 + 2)(A;2 -q2 + 1)(P2 + g2 - A:2) 
^ (A;2-Z2 + l)(A;2-/2 + 2) 
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{k-l~ki + li + k2-l2 + 2) 



{k-l + ki-li + k2-l2 + ^y 




(/2 + l)(g2-/2)(P2 + g2-/2 + l) 
(fc2-/2)(A;2-/2 + l) 



X 



{ki-li + k2-l2-k + I) 



-l + ki-h-k2 + l2 + 2) 



and 



(A; + 2) 



(31) 



^ {k-q + l){p + q-k)' 



The SU{3) CGC are found by using Equation H. The corresponding SU{2) 
coefficients are determined through Equation |1^. 

3.2 Algorithm for Isoscalar Factors 

The use of the above recursion relations to generate isoscalar factors is far from 
trivial, and a detailed description of the algorithm is given in |2^. A summary 
of that procedure is given here. 

The ISFs for coupling to a state of highest weight can be thought of as 
occupying lattice sites in an irregularly-shaped volume in a three-dimensional 
space. With the restriction provided by hypercharge conservation, these ISFs 
can be thought of as functions of s = Ji + /2, fci, and ^2. The two recursion 
relations above for the shw states each have two terms with the same s value, and 
two others with s values smaller by one s—1. All ISFs with the maximum allowed 
s value for a particular coupling are related, therefore, by a two-term recursion 
relation. Setting one of these equal to one (to be fixed via normalization later), 
all others are simply constructed. Because of the shaped of the lattice, one 
can in most cases find a particular single ISF with s = Smax — 1 which can 
be determined by those with s = Sma,x'- once it is determined, all others with 
the same s value can be evaluated through the complete four-term recursion 
relations based on values already known. Repeating this logic, one moves to a 
single ISF at the next lowest s value from which all others with the same s can 
be determined, continuing until all ISFs are known. Normalization requires the 
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sum of the squares of all the ISFs so determined must be 1, so each is multiplied 
by the appropriate normalization factor. 

When dealing with ISFs of degeneracy greater than one, it is necessary to 
make several choices of ISF values in order for the algorithm to succeed. Which 
are chosen, and what values they are given, determines the choice of degeneracy 
resolution. The criteria mentioned in an earlier section is equivalent to choosing 
ISFs for the d — 1 highest values of s to be zero, and following the algorithm of 
the previous paragraph beginning with the s = Smax — d + 1 values. A second set 
begins with the s = Sma,x — d+2values, and the further condition of orthogonality 
of these ISFs to the previous ones. Subsequent sets of ISFs are formed likewise 
by successively starting with higher values of s, and enforcing the condition of 
orthogonality to those previously determined. 

Regardless of degeneracy, there are particular combinations of representa- 
tions whose ISFs cannot be determined as above due to a failure in the ability 
to step to lower values of s. In every such case, the algorithm succeeds for the 



conjugated ISFs, and the exchange symmetry relation of Equation |T7| can be 
used to find the original values. 



4 Program Structure 

We define a variable type fraction as an array of two long integers, n and m, 
to represent a number of the form ^n/m. All SU{3) CGCs and ISFs can be 
exactly represented this way. The routines that manipulate variables of the type 
fraction are 



reducef racO 
setequalf racO 
addf racO 
subf racO 
multf racO 
imultf rac 
divf racO 



Reduces a fraction to its lowest form. 

Assigns one fraction to be equal to another. 

Adds two fractions. 

Subtracts two fractions. 

Multiplies two fractions. 

Multiplies a fraction by a long integer. 

Divides two fractions. 
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absf racO 

f actorialf rac 
equalf racO 
comparef racO 

addirrat () 



Changes a fraction to its absolute value and returns 
its sign. 

Finds the factorial of a fraction, if it is integral. 
Returns true if and only if two fractions are equal. 

Returns true if and only if the first fraction is greater 
than the second. 



Adds two numbers of the form ^n/m. Since the co- 
efficients of SU (3) always take this form, no incom- 
patibilities are expected. An error message is written 
to the log file if it is not possible to find 



M ~ \ nil \ m2 



(32) 



Global variables are 



pi, ql, p2, q2 

numprod 

p[], q[] 
names [] 
texnames [] 



The identities of the factor representations 

and (P2,g2)- 

The number of product representations in the 
Clebsch-Gordan series. 

The identities of the product representations {pi,qi). 

The names of the representations in the product. 

The BTgX names of the representations in the 
product. 

ISFtable [] [] [] [] [] [] Array which holds the isoscalar factors. Its elements 

are of the type fraction. 

SHWtable [] [] [] [] Array which holds the isoscalar factors for the states 

of highest weight. Its elements are of the type 
fraction. 

step_down [] [] Array which holds certain ratios of ISFs for doshwO . 
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flag 



errcount 



An integer flag used by dosliceO to control 
doshwO. 



Cumulative count of errors in the program's 
execution. 

nameconv, labelconv Integers denoting the user's choices for labelling con- 
ventions (sec below). 

maxr, minkl, maxkl. Limits used in dimensioning the various arrays, 
maxll, mink2, 
maxk2 , maxl2 

Other subordinate routines included are 



sizerepO 

ipower () 
factorial 
gcdO 

isqrtO 
yiiO 

digitstringO 
numstringO 
namerepO 
shwO 

initlims 



Actually a macro which returns the dimension of a 
representation. 

Raises a long integer to a power. 

Returns the factorial of a long integer. 

Returns the greatest common devisor of two long in- 
tegers. It is employed by reducef rac(). 

Approximates the square root of a long integer with 
another long integer. 

Gives the hypercharge and isospin of a state with 
given k, I, m. 

Gives a character representing a one-digit integer. 
Used by numstring. 

Returns a character string that represents an integer. 
Used by namerepO. 

Gives the text name and BTgX name of a represen- 
tation. Naming conventions are explained below. 

Gives the hypercharge and isospin of the highest- 
weight state of a representation. 

Sets the values of minkl, maxkl, maxll, mink2, 
maj?;k2, ma:xl2, and maxr. 



14 



inittableO 

kill_table() 
put entry 



get entry 

initshwO 

killshwO 
putentrySHWO 
getentrySHWO 
printf ractionO 

Printline 
numstatesi 

numstatescO 

getconvs 



Allocates memory space for ISFtable and initializes 
its values to zero. 

Releases the memory used by ISFtable. 

Places an entry into ISFtable. Because of liy- 
percharge conservation, one index is redundant 
for ISFtable; this is used by putentryO and 
getentryO to reduce the size of the array for 
ISFtable. 

Checks to see if indices defining an ISF are all in 
range, and if so, returns an entry from ISFtable. 

Allocates memory space for SHWtable and initializes 
its values to zero. 

Releases the memory used by SHWtable. 
Places an entry into SHWtable. 
Returns an entry from SHWtable. 

Writes a fraction in IXTgX format to one of the out- 
put files. 

Writes a line to both cle.tex and iso.tex. 

Counts the isomultiplets of given hypercharge and 
isospin in the product. 

Counts the states of given hypercharge and third 
component of isospin in the product. 

Prompts for and inputs the naming and labelling 
conventions (see below) from the user. 



The central routines are 



su2clebsch() 



degeneracy 



Calculates an SU (2) Clebsch-Gordan coefficient with 
Equation [1^. 

Returns the number of times a representation ap- 
pears in the Clebsch-Gordan series. 
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series 

texphasesO 

normhor 

normvert 

checktableO 
checksymO 



textableO 



coef xO 



llstepO 
klstepO 
doshwO 

dosliceO 



Finds the representations in the Clebsch-Gordan se- 
ries. Also sets the value of numprod. 

Tabulates the symmetry phases and 1^3 into 
iso . tex. 

Checks that each row (for each isospin) of the ISF 
table has unit norm. 

Checks that each column of the ISF table has unit 
norm. 

Calls normhor and normvert () . 

Checks that the shw ISFs are symmetric under the 
exchange of ri and r2, in the cases where pi = p2, 
ql = q2. It compares only the absolute values of the 
ISFs. It is called by dosliceO. 

Loops over the quantum numbers of the factors 
and products and tabulates the isoscalar factors 
into iso. tex. Calls su2clebsch and tabulates the 
Clebsch-Gordan coefficients into cle.tex. 

Returns a coefficient of the recursion relations in 
Equations |21|-|3T|. Here x is one of al, a2, a3, a4, 



bl, b2, b3, b4, cl, c2, c3, c4, dl, d2, dS, alpha, 
beta. 

Steps up or down in li for doshwO. 

Steps up or down in ki for doshwO. 

Fills the parts of ISFtable for the highest-weight 
states of the product representations. 

Fills ISFtable for a given {pi, qi) in the product. 
Uses the algorithm given in the previous section. 
Calls doshwO to handle the states of highest weight. 
Calls checksymO for cases of pi = p2, ql = q2. 
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mainO Prompts for pi, gi, p2, q2- Calls getconvsO 

to input labelling and naming conventions. Uses 
degeneracy to determine the identities of the prod- 
uct representations and calls dosliceO for each. 
Calls textableO to generate output. Also calls 
normhorO and normvert () . 



5 Using the Program 
5.1 User Input 

The user, upon initiating the code, is queried to supply the identities of the 
factor representations in the form (pi, gi), (p2, Q'2), where each is an integer. The 
naming and labelling conventions must also be chosen. The choices for naming 
convention are 



1. Representations are named as in |2T|, which is popular among particle 



physicists. Representations are named by their dimension and bars are 
added based on triality 

t = {p — q) mod 3. (33) 

Representations with t = 1 are unbarred, and those with t = 2 are barred. 
One exception is that 6 = (2, 0). Representations with t = have bars if 
and only if g > p. 

2. Representations are named by their dimension and bars are added if and 
only if g > p. 

3. Representations are named by (p, g). 

In choices 1 and 2 above, there is an ambiguity when two representations with 
different (p, g) have the same dimension. We distinguish them using primes. In 
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the naming scheme 1 the lowest-lying representations with same dimensions are 



15 = (2,1), 

105 = (6,2), 

120 = (3,5), 

195 = (9,2), 

210 = (4,6), 

231 = (2,10), 

The choices for labelling are 



15' = 


(4,0), 


105' = 


(13,0), 


120' = 


(1,9), 


195' = 


(1,12), 


210' = 


(19,0), 


231' = 


(0,20). 



120" = (0,14), 



1. States are labelled by hypercharge and isospin y, i, i^. Isomultiplets are 
labelled by y and i. 

2. States are labelled by projection quantum numbers k, I, m. Isomultiplets 
are labelled by k and I. 

5.2 Error Messages 

The program generates three files during its execution. The log file (logf ile 
= su3 . log) contains messages and reports the progress of the program. See 
Sample 1 below. Any error messages are written to the logfile. The most 
important errors are also written to the console. They are: 

WARNING: insufficient memory 

indicates that the free memory in the computer is too small to allow dimension- 
ing of the needed arrays. 

WARNING: integer overflow 

indicates that the integer mathematics has generated numbers larger than the 
"long" integers defined by the computer. 

WARNING: Table is not horizontally normal 
indicates that the table of ISFs has nonnormal rows. A row is given by 
kl, 11, k2, 12, and the isospin of the states in the product. 

WARNING: Table is not vertically normal 

indicates that the table of ISFs has nonnormal columns. A column is given by 
R,k,l. 
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In addition, for cases in which pi — p2 and ql — q2, a, quick check is made 
of the exchange symmetry. If a nonsymmetric entry is found, the message 
WARNING: n=4 SHW ISFs are not symmetric 

(for example) appears. 

These errors indicate that a serious problem has occurred. Other error 
messages are written only to logf ile and are self-explanatory. 

5.3 Output 

The ISFs are written in BTgX format to the file isofile = iso.tex. In it, 
representations are named according to the naming choice made at the beginning 
of execution, and states are labelled according to the labelling choice. Multiply 
degenerate representations are distinguished by subscripts. A square root is 
assumed to appear over the unsigned part of each entry. Thus, 





R 


ri r2 


Y 




I 


yi ii Z/2 i2 


±C 



(35) 



means that the isoscalar factor 

F(R, Y, I : ri, yi, h; r2, y2, ^2) = ±VC (36) 
in the {y, i) labelling convention. In the (/c, /) convention, this would appear as 





R 


ri ® r2 


k 




I 


ki h k2 h 


±C 



(37) 



The symmetry factors ,^1 and .^2 are also given in this file. For example, 



ri ® V2 


R 




+ 


6 
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(39) 



means that 

ei(R:ri;r2) = +1, 
e2(R:ri;r2) = -1. 

The Clebsch-Gordan coefficients are written in IMgX format to the ffie 
defile = cle.tex. Again, a square root is assumed over the unsigned part of 
each entry. In the {y,i,iz) labelhng convention, the format is 



(40) 





R 




Y 


ri (g) r2 


I 




Iz 


Vl H Hz 1/2 i2 i2z 


±C 



for 

(R r / /^|ri yi ii ii^ ra 2/2 ^2 ^22) = ±VC. 
In the {k,l,m) notation this is 





R 




k 


ri (g) r2 


I 




m 


ki h nil k2 h ^2 


±C 



(41) 



(42) 



Samples 2 and 3 below show examples of iso.tex and cle.tex after pro- 
cessing by lATgX. 



6 Sample Output 

Presented here are sample output ffies from a successful run. The product 
computed is 15 ® 8 = (2, 1) ® (1, 1), which has a doubly degenerate 15. 
In this run, we have chosen the naming convention of and the labelling 
convention of hypercharge and isospin {y, i). Sample 1 is the log file su3.log. 
It echoes our naming and labelling conventions, informs us of the program's 
progress, and tells us that the table of isoscalar factors is properly normalized. 
Sample 2 is the IMgX ffie containing the isoscalar factors. The Clebsch-Gordan 
coefficients are in the IMgX file of Sample 3. Notation for the tables is explained 
in the previous section. 
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SAMPLE CAPTIONS 



Sample 1: Log file (su3 . log) from a sample run calculating 15 8 = 
(2,1)0(1,1).. 

Sample 2: Output file (iso .tex) from a run calculating 15 8. These 
are the isoscalar factors. Notation is explained in Section 4. 
The output has already been processed by IMgX. Only the 
first page is shown, in order to save space. 

Sample 3: Output file (cle .tex) from a run calculating 15 8. These 
are the Clebsch-Gordan coefficients. Notation is explained 
in Section 4. The output has already been processed by 
l^TgX. Only the first page of the file is shown, in order to 
save space. 
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SU3 : beginning 



getconvs: naming convention is of Slansky 
getconvs: labelling convention is (y, i, i3) 

SU3: calculating 15 x 8 



sua 


working 


on 


3 with degeneracy = 1 


SU3 


working 


on 


6* with degeneracy = 1 


SU3 


working 


on 


15' with degeneracy = 1 


sua 


working 


on 


15 with degeneracy = 2 


SU3 


working 


on 


24 with degeneracy = 1 


SU3 


working 


on 


42 with degeneracy = 1 



checktable: table is horizontally normal 
checktable: table is vertically normal 

SU3: writing output files 

SU3 : ending 
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Isoscalar Factors for 15 (8) 8 



15 (8) 8 
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15' 15i 


152 


24 


42 
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+ 




+ 
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+ 




+ - 


+ 
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24 




42 
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4 1 1 1 


-1 


4 1 1 1 
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3 2 




3 2 







6 24 
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1 1 1 1 
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5 "5 
4 1 
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42 
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400 
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16 
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27 














1 
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529 
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Clebsch-Gordan Coefficients for 15 (8) 8 
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